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In multiplex networks with a large number of layers, the nodes can have different activities, 
indicating the total number of layers in which the nodes are present. Here we model multiplex 
networks with heterogeneous activity of the nodes and we study their robustness properties. We 
introduce a percolation model where nodes need to belong to the giant component only on the 
layers where they are active {i.e. their degree on that layer is larger than zero). We show that 
when there are enough nodes active only in one layer, the multiplex becomes more resilient and the 
transition becomes continuous. We find that multiplex networks with a power-law distribution of 
node activities are more fragile if the distribution of activity is broader. We also show that while 
positive correlations between node activity and degree can enhance the robustness of the system, 
the phase transition may become discontinuous, making the system highly unpredictable. 

PACS numbers: 


I. INTRODUCTION 

Multilayer networks m are formed by several inter¬ 
acting networks. They describe a large variety of com¬ 
plex systems. Examples are found in social [1], tech¬ 
nological, communication, transportation systems mi 
but also in biological networks of the cell, or in the brain 
[9Hn]. In the last fifteen years a lot of attention has 
been devoted to understand the interplay between struc¬ 
ture and dynamics in single networks In recent 

years, it has become clear that most of the networks are 
not isolated and that for understanding the function of 
complex systems as different as complex infrastructures 
or the brain it is important to investigate the role of their 
multilayer structures [HE]. In particular, ample debate 
has been devoted to characterize the robustness of multi¬ 
layer networks dniElHSO]. It has been shown that, in the 
presence of interdependencies, the robustness of a mul¬ 
tilayer network can be significantly affected. Multilayer 
networks can be much more vulnerable to random dam¬ 
age with respect to considering only their single layers 
taken in isolation dlHIll. 

In the presence of interdependencies, the notion of mu¬ 
tually connected component has been introduced, mean¬ 
ing that each pair of nodes in the mutually connected 
component must be connected by a path on each and ev¬ 
ery layer, internal to this component. This definition is 
motivated by the fact that in these interdependent sys¬ 
tems a node is not functional if any of its interdependent 
nodes in the other layers is not functional. Therefore, the 
largest (giant) mutually connected component (MCGC) 
describes the robustness of the system. This component 
has a discontinuous phase transition as a function of the 
initial random damage inflicted to the nodes of the net¬ 
work, and close to this transition the system is affected 
by dramatic cascades of failure events. This transition 
can be studied on multilayer networks of different nature, 
including multiplex networks [IMZl and network of net¬ 
works EDHUj. Network of networks are formed by dif¬ 


ferent interacting networks, as the molecular networks in 
the cell, in which every node is a different type of biolog¬ 
ical molecule. Multiplex networks [g H EIHSe], instead, 
are multilayer networks in which the same set of nodes 
interact through different networks. Example of multi¬ 
plex networks are social networks, in which individual are 
connected by different types of social ties, transportation 
networks, like airport networks in which each airport is 
connected to other airports though connections of differ¬ 
ent airline companies, or collaboration networks in which 
scientists collaborate on different topics and eventually 
cite each other. Several works have focused on modeling 
multilayer networks [32H36j . In particular a very useful 
approach employs statistical ensembles which are able to 
generate a large variety of multiplex network topologies 
with desired level of structural correlation [311 EH 155] . 
In networks of networks as well as in multiplex networks 
the possible correlations existent in this structure can 
strongly affect their robustness properties m Eg EH EH). 

An interesting result has recently shown that multi¬ 
plex networks are characterized by the fact that not all 
the nodes are connected in every layer. In fact, many 
networks have been shown [H to have heterogeneous ac¬ 
tivity of the nodes. The activity of a node is given by 
the number of layers in which the node is at least con¬ 
nected to another node. The activity of the nodes has 
been found to be broadly distributed in a variety of mul¬ 
tiplex network data sets (H- Moreover, the activity of 
the nodes has been seen to correlate in average with the 
degree of single node in given layers, meaning that on 
average nodes that are present in many layers have also 
typically high degree in the single layers in which they are 
active (H • The heterogeneous distribution of activities of 
the nodes and its relation with the mean degree in sin¬ 
gle layers is frequently observed in real-world multiplex 
networks and may encode relevant information |H. Be¬ 
sides, multiplex networks with heterogeneous activities 
can be useful theoretical tools to address other problems 
like modular structures in single networks [5|. 


2 


Here we characterize multiplex networks ensembles 
with heterogeneous activities of the nodes and we investi¬ 
gate their robustness properties, assuming mutual inter¬ 
dependencies among each layer. This ensemble of mul¬ 
tiplex networks with heterogeneous activity of the nodes 
can be used to model realistic multiplex network struc¬ 
tures, on top of which dynamical processes may occur. 
Moreover, ensemble of multiplex networks with heteroge¬ 
neous activities of the nodes can be used to estimate the 
role that correlation have on their robustness properties. 
Here we find that heterogeneous activity of the nodes 
can decrease the robustness of networks in the presence 
of interdependencies. Nevertheless the correlation be¬ 
tween the activity of the nodes and their degree within 
single layers has the opposite effect and can improve the 
robustness of multiplex networks. 

II. THE MULTIPLEX NETWORK ENSEMBLE 
WITH HETEROGENEOUS ACTIVITY OF THE 
NODES 

A multiplex network of N nodes i = 1,2,... ,N and 
M layers a = 1,2,... ,M is completely specified when 
the M adjacency matrices of elements a“- = 1 if node 
i is linked to node j in layer a, and otherwise a“- = 0 
are given. Every node of the multiplex network is la¬ 
belled as {i,a), indicating that is the i-th. node in layer 
a. The replica nodes of the node {i, a) are defined as all 
the nodes labelled as {i,a') in layers a' ^ a [ 53 ]. Inter¬ 
estingly, it has been observed from data |7] , that in many 
networks not all the nodes are active {i.e. are connected 
to at least another node) in each layer. Let us define the 
activity Bi of a node i as the number of layers where 
node i has non-zero degree. The activities of the nodes 
are broadly distributed |7] and they can be fitted by a 
power-law P{B) ~ B~^ with S € [1.5, 3.0]. This implies 
that for some multiplex networks the bipartite network 
between nodes and layers described by the activity adja¬ 
cency matrix can be either dense S < 2 or sparse 5 > 2 
but the typical number of layers in which a node is ac¬ 
tive is always subject to unbound fluctuations. In order 
to characterize fully the activities of the nodes in each 
layer, from the M adjacency matrices a“ it is possible 
to construct a N x M activity matrix b of elements bi a 
(Fig. §. This matrix can be viewed as an adjacency 
matrix between nodes and layers indicating if node i is 
active in layer a {bi^a = 1 ) or not {bi^a = 0). The activity 
Bi of a node i can be therefore expressed in terms of the 
matrix b as in the following, 

M 

Bi = '^bi^a- ( 1 ) 

a—1 

The layer activity Na has been defined in mm and is 
given by the number of nodes active in layer a, i.e. 

N 

A^ck 'y ( bi^a- 

i=l 




FIG. 1: Sketch of a multilayer network with node activity: 
nodes are not active in every layer (left), so the system can 
be represented by a bipartite graph associating each node to 
the layers where it appears (right). 

It is therefore possible to construct ensemble of multi¬ 
plex networks with heterogeneous activity of the nodes 
in the subsequent manner: first we construct the network 
between layers and nodes described by the adjacency ma¬ 
trix bi a indicating if node i is active in layer a, then we 
can construct in each layer a network between the active 
nodes of the layer with a degree distribution P°^{k). In 
the case in which P°‘{k = 0) = 0, we have that a node is 
active in a given layer if it is connected at least to another 
node in the same layer, i.e. 

bi a = 1 , a“ : (3) 

where indicates the Kronecker delta. 

In order to construct the activity matrix bi^a we can 
consider a microcanonical ensemble in which the activi¬ 
ties of the nodes Bi and the layer activities Na are fixed. 
This ensemble, called also the configuration model of a 
bipartite network, can be constructed by maximizing the 
entropy S of the activity networks given by 

5 = -^P(b)logP(b) (4) 

b 

where H(b) is the probability of a given activity network 
in the ensemble under the constraints that the node activ¬ 
ity and layer activity are kept constant. In this ensemble 
the probability of a matrix b is given by 

N / M \ M / N \ 

^(b) = ^nM ’ (5) 

i—1 \q(— 1 / ct—1 \i—1 / 

where Z is the normalization constant also called the 
partition function of the ensemble. In this ensemble the 
probability pi^a that a node i is active in layer a is ex¬ 
pressed in terms of the Lagrangian multipliers \i, Wq, 
ie. 

g-Ai-Wc 

The Lagrangian multipliers are fixed by the conditions 
^ ^ Pi,a 

i 

y '. Pi,a — Pi 


( 2 ) 


( 7 ) 
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Finally the entropy |55HiD] S of the ensemble of activity 
networks is given by 

S = S-il (8) 


III. MUTUALLY CONNECTED COMPONENT 
IN A MULTIPLEX NETWORK WITH GIVEN 
DISTRIBUTION OF ACTIVITIES OF THE 
NODES 


where S and are given by 


5 = 

n = 


- X! + (1 

ict 





-El- 


NJ 


N, 




We assume here that bipartite networks between nodes 
and layers, characterized by the adjacency matrix bi^a 
is uncorrelated, i.e. we assume that the activities Bi 
of the nodes are not correlated with the sizes of active 
nodes on the layers where the nodes are active. In 
this hypothesis the probabilities pi^a are proportional to 
the product of Bi and Na which are the degrees of the 
mentioned bipartite network, and the following condition 
on the the maximal activity of the nodes B^^x and the 
maximal layer activity Nmax must be satisfied, i.e. 


Bmax^max 


< 1 . 


(9) 


This condition is the condition that in a bipartite network 
corresponds to the one imposing a structural cutoff on 
single uncorrelated networks (networks without degree- 
degree correlations) [33. In this case we have the simple 
expression 


Pi,oc 




( 10 ) 


with Np = E^Bi = {N)M = {B)N. 

As we said, this argument regards lack of correlation 
between activities and number of active nodes on the 
layer. In Sec. IV, we will examine a different type of 
correlation, namely the one between activities and the 
degree distribution on the layer. 

Once the activity network is constructed, in order to 
construct the multiplex network, we assign to each ac¬ 
tive node of layer a a degree within the layer. If the 
activity of the nodes are uncorrelated with the degree in 
each layer the degrees kt^a of the nodes i in layer a are 
drawn form the degree distribution P°‘{k) and the net¬ 
works of each layer a are generated by the configuration 
model with degree distribution P°^{k). In other words, 
the probability of the multiplex network degree sequences 
k = {ki^a}i=i,2...,N-a=i,2....,M IS given by 


We consider the mutually connected giant component 
(MCGC) in a multiplex network with given distribution 
of activities of the nodes, described by the ensemble of 
multiplex networks introduced in section [II] 

The layers are interdependent, meaning that each node 
active in a given layer a is interdependent on its replica 
nodes in all the layers /3 where the node is active. In 
particular we will assume that a node (i, a) active in layer 
a belongs to the mutually connected giant component if 

(a) at least one neighbor (j, a) of node (f, a) belongs 
to the mutually connected giant component; 

(b) in each layer j3 ^ a where the node i is active, 
at least one neighbor (j, (3) belongs to the mutually 
connected giant component. 

On a locally tree-like multiplex network, this can be eas¬ 
ily encoded in a message passing algorithm determining 
if a node (i, a) belongs to the mutually connected gi¬ 
ant component [D ng i3oi l4T] . The indicator function 
Sia indicates if a nodes (i, a) belongs (Sia = 1) or not 
(Sia = 0) to the mutually connected giant component. 
This indicator is determined by a set of “messages” 
that each node (i, a) active in a layer a sends to the 
neighboring nodes (j, a) in the same layer. Each mes¬ 
sage indicates if the node (f, a) belongs = 1) 

or not (o'jU.j = 0) to the mutually connected component 
when the link to the node (j, a) is removed. 

Here our goal is to characterize the size of the mutually 
connected giant component as a function of the probabil¬ 
ity 1 — p that random nodes are damaged in the network. 
In order to characterize the damage initially inflicted to 
the network we indicate with Sj = 0,1 if a node has been 
damaged (si = 0) or not {si = 1) in the multiplex net¬ 
work. Given the definition of the mutually connected 
component, the message from a node (i, a) to a 

node (j, a) both active in layer a is therefore given by 




1 - n 


xn 

/3 

/3/a 


i-k/3 n 

ledfiii) 


( 12 ) 


P({k}|b) = n [P“(fc“)&.,„ + Sikf,0)il - 6,,„)]. (11) 

i,a. 

Instead, if there is a correlation between the degree of the 
nodes within a layer and the activity of the nodes, then 
we need to draw the degree of the nodes in each layer 
from a probability Pg^i^ki) which is a function of their 
activity Bi. 


where da{i) indicates the nodes i that are neighbors of 
node i in layer a and the expression da(i)\j indicates all 
the set of all the nodes belonging to da{i) except node j. 
Therefore the message crjkyj is equal to one, if the node 
has not been damaged (si = 1), if at least a message 
coming from a neighbor node {£',a) different from (j, a) 
is positive (first factor in the square brackets), and if in 
all the layers /3 where node i is active = 1), there 
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is at least one neighbor (€,/3) of node (i,/?) sending a 
positive message to node (i, j3) (product over /3). Finally, 
the indicator function that node i active in layer a, is in 
the MCGC is given by 


— 


1 - n (1 




n 


/3/a 


1 “ Kp n (1 - ct: 






(13) 


Vi,Cl — 


B^Nq, 

{B)N 


< 1 , 


we have 


1 < < 


{B)N 
B, ’ 


(14) 


(15) 


{B)N 

M 


» 1 . 


(16) 


This means that, for example, for finite (B) the number 
of layers is much smaller than the total number of nodes. 

In order to characterize how the size of the mutually 
connected giant component depends on the distribution 
of activities P{B), here we consider the ensemble of mul¬ 
tiplex networks with given activities of the nodes de¬ 
scribed in section HH where we assume additionally that 
the damage occur on a node i with probability 1 — p, i.e. 


N 




(17) 


i=l 


Let us observe now that the probability that a node i 
active in layer a has activity Bi = B is given by 

P{Bi = B,bi^a = 1 ) 


P{B, = B\h,a. = 1) = 


= P{B) 


P{h,o. = 1) 

BN, 


{B)N 
P{B)B 




BN„ 

{B)N 


where we have assumed that pi^a is given by Eq. (lOl. 
Assuming that all the layers have the same activity 
Na = N MjS, and the same degree distribution P{k), 
by averaging the messages given by Eq. (12) over the en- 


This indicator function equals one, if the node (i, a) has 
not been damaged (s^ = 1), if at least a message coming 
from a neighbor node {P,a) of layer a is positive, and 
if in all the layers /3 where node % is active = 1), 
there is at least one neighbor (£, /3) of node (I, /3) sending 
a positive message to node (i, j3). In this section we only 
consider uncorrelated activity networks in order to allow 
for a theoretical treatment of the percolation properties 
of the multiplex networks with heterogeneous activity of 
the nodes. Moreover we assume that the number of nodes 
Na is large in every network a of the multiplex network. 
Given that from Eq. (10) we have 


semble of multiplex networks with given activities of the 
nodes, we get the probability cr“ = that by fol¬ 

lowing a link in layer a we reach a node in the mutually 
connected component, obtaining cr“ = tr and 

^ = pE^^^[1-Gi(1-^)H1-Go(1-<7)]^-^(19) 

where the generating functions Go{x) and Gi(x) are 
given by 


Go(a:) 

Gi{x) 






.k-1 


( 20 ) 


Which can be derived as in the following. First we 
define cr“ as 

B 

where (...) indicates the average over the network, in the 
large network limit. Using Eq. (HI we have 

a = p[l- Gi(l - a)] Y, P{B^ = B\b,^^ = 1) 


where the condition 3> 1 is required to study the 
percolation transition of the MCGC. Since the value of 
the activity of the nodes Bi < M we find the condition 



ledj3{i) 


( 21 ) 


= p[l _ Gi(l - a)] ^ - Go(l - a)f-\ 


getting Eq. (191 in the tree-like local approximation. 


Using a similar derivation, under similar assumptions, 
the fraction S'" = (Sia) = S of the N nodes that are in 
the mutually connected component in layer a is given by 




B 


Hgp[i _ Go(l - a)] [1 - Go(l - <y)f-H22) 


These equations can be studied in detail as a function 
of the degree distribution in each layer, and the node 
activities distribution P{B). In particular, Eq. (19) can 
be expressed as. 


a = p[l- Gi(I - a)]K (1 - Go(l - a )), 


(23) 


where we have indicated by K{x) the following generat¬ 
ing function 


(18) 






(24) 
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Similarly 3°“ = S can be written as 

S=p[l- Go(l - a)]K (1 - Go(l - a)). 


(25) 


Let us consider now Poisson networks identical on each 
layer: P{k) = e~‘^c^(k\ .Then we have Go{z) = Gi{z) = 
g-c(i-z) equation for cr becomes h{x) = 0 where 


h{x) = X — p{l — e ^)K{l — e ^), 


(26) 


and X = ca, p = cp. A discontinuous transition occurs 
only if h(x) has maximum, so we need to impose the 
condition h{x) = h'{x) = 0, where 

h'{x) = 1 - pe-^ [K {1 - e-^) + {1 - e-^)K'{1 - e-^)] 

where the function K'(x) is given by 


K'{x) = Y. 


(B) 


-P{B)x 


B-2 


(27) 


Using this equation to isolate p, and substituting into 
h{x) = 0, we get for a; ^ 0 


A:'(l-e-“) _ 1 

AT (1 — e~^) X 1 — e~ 


(28) 


As a general remark, we expect in this model that the 
nature of the percolation phase transition will be depen¬ 
dent on the fraction of replica nodes with activity B — 

In fact these nodes do not have any interdependency on 
replica nodes in other layers. Therefore a high density of 
nodes with activity B = 1 favors the continuous phase 
transition. In the extreme case in which all the nodes 


have activity B = \ Eq. (261 is in fact reduced to the 


equation determining the giant component of a single net¬ 
work. In fact, the generating function K(x) defined in 
the last equation of ( p4| ) is equal to one if all the nodes 
have activity B = 1. When the not all the nodes have 
activity B = 1, the non trivial functional form of K{x) 
is determined by the nodes that have activity B > \. 
Therefore we expect that the density of such nodes is 
sufficiently high the transition is continuous while other¬ 
wise, we expect a discontinuous phase transition. This 
expected phenomenon can be related with a similar ef¬ 
fect generated by the partial interdependence in multi¬ 
plex networks where all the nodes are active in any layers 

[niiiT]. 

In the following we provide evidence for this observa¬ 
tions by considering two different distributions of activi¬ 
ties P{B): a scale-free activity distribution and a Poisson 
activity distribution. 


A. Scale-free P{B) distribution 


Let us consider a power law activity distribution 
P{B) = B -^with B > Bmin- In this 


case, the function h{x) is given by Eq. ( [26| where 

E OO 

B^Brr 


K{x) = 


B1-SxB-1 


E OO 

B^l 


B 


l-<5 


(29) 



FIG. 2: Phase diagram of a multiplex network formed by 
Poisson layers in the case of a power law activity distribu¬ 
tion with exponent 5. The solid line indicates the continu¬ 
ous phase transition obtained for Bmin = 1, and the dashed 
line indicates the discontinuous phase transition obtained for 
Bmin 2. 

10 r- 
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FIG. 3: Jump in Xc = ca^ = cSc, where Sc is the size of 
the MGGG at the discontinuous transition for a multiplex 
network formed by Poisson layers with mean degree c, in the 
case of a power law activity distribution with exponent <5 and 
Bmin — 2 . 



The position of the discontinuous transition can be cal¬ 
culated by solving h{x*) = h'{x*) = 0, with x* > 0. 
The discontinuity of the transition is due to the finite 
value of the solution at x* > 0, as one can see that 
S > 0 from Eq. (251, where Go{z) is a simple exponen¬ 
tial. The position of the continuous transition can be 
instead calculated by solving h{0) = h'{0) = 0. While 
the tricritical point, if it exist, can be found by setting 
h{0) = h'{0) = h''{0) = 0. (Note, that given the func¬ 
tional form of h{x), given by Eq. (261, the condition 
h{0) = 0 is always satisfied.) 

Fig. |2] illustrates the phase diagram calculated for 
Bmin = 2 and Bmin = 1- For Bmin = 1 we have a contin¬ 
uous phase transition for S > 2. For Bmin = 2 we find a 
discontinuous phase transition line across all the values of 
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d > 2. The jump in the size of the MCGC Xc = cac = cSc 
at the discontinuous transition for Bmin = 2 (shown in 
Fig.§ diverges as <5 approaches the value 2 (while at the 
continuous transition there is no jump). Both transition 
lines for Bmin = 1 and Bmin = 2 diverge as S approaches 
the value 2, indicating that as <5 decreases the multiplex 
network becomes more fragile. As 6 decreases, the frac¬ 
tion of nodes with large activity (activity hubs) becomes 
important. These nodes, to be in the mutually connected 
giant component they must be connected in each layer, 
and if they are damaged they affect the connectivity of 
multiple layers. Therefore they are at the same time more 
fragile than nodes with smaller activity and are able to 
affect the multiplex connectivity across more layers. As 
a result, as 6 decreases, random damage affects more and 
more layers and the multiplex network becomes less ro¬ 
bust. 


B. Poisson P{B) distribution 

For activities that are Poisson distributed P{B) = 
fi^e~^/B\, we have 

h{x) = X — p(l — e~^) exp [—^(e”"^)] . (30) 

By setting h{x*) = h'{x*) — 0, we can find a line 
of discontinuous transitions that ends at a critical point 
defined by h{0) = h'{0) = /i"(0) = 0 characterized by 

p = Pt = Ve 1M8, p = pj’ = ^= 0.5 (31) 

In the multiplex network, when p < px the percolation 
transition is continuous. Instead for p > px the percola¬ 
tion transition is discontinuous. In Figure]^ we show the 
phase diagram. At p = 0, the continuous phase transi¬ 
tion is at p = pc = 1 as expected for independent layers. 
As the level of interdependence becomes increasingly sig¬ 
nificant, the continuous transition becomes discontinuous 
(for p > px). This can also be seen from the size of the 
discontinuous jump Xc = cuc = cSc, that approaches zero 
at the tricritical point px and is shown in Figure 
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FIG. 4: Phase diagram of a multiplex network formed by 
Poisson layers with identical mean degree c and with a Poisson 
activity distribution with mean p. The blue solid line is a line 
of continuous phase transitions between a percolating and a 
non-percolating phase, the red dashed line represents instead 
the line of discontinuous phase transitions. 
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0l - ^^^^^- 

0 1 2 3 4 5 

FIG. 5: The jump in Xc = cSc where Sc is the size of 
the mutually connected giant component at the discontinu¬ 
ous transition point, is shown for a multiplex network formed 
by Poisson layers with identical mean degree c and with a 
Poisson activity distribution with mean p. 


IV. CORRELATIONS BETWEEN THE 
ACTIVITIES AND THE DEGREES OF THE 
NODES 

In this section, we consider the role of correlations be¬ 
tween the activities of the nodes and their degrees in the 
layer in which they are active. The higher is the activity 
of a node, the more fragile is the node, but the higher is 
the degree of the node within the layer the more robust it 
is. Therefore, here we want to characterize the effect that 
the correlations between the activity and the degrees of 
each node have on the robustness properties of the entire 
multiplex network. We start from the message passing 
Eqs. (12) and (13). We assign to each node i an activity 


Bi from a P{B) distribution. Then, we assume that the 
matrix of activities is given by b with probability P{h) 
given by Eq. (1^ with pi ^ given by Eq. (10). Finally, we 
need to extend Eq. 0 for the degree sequence to the 
correlated case. As node degrees on each layer are cor¬ 
related with node activities, the probability P({k}|b,B) 
that a multiplex network has degree sequences {k}i a, 
given the activity matrix b and the activity sequence 
{Bi}i=i^2...,N, is 


P({k}|b, B) = [| [Pb,( knbi.„ + 5(kr, 0)(1 - bi,„)] 


( 32 ) 
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FIG. 6: Phase diagram of a multiplex network in the case 
of an activity distribution P{B) oc B~^ for B £ [1,100] and 
degree within each layer correlated with the activity of the 
node according to Eqs. (331 — (34l. From top to bottom, 
the dashed lines indicate discontinuous phase transitions for 
a = 0.2 and for a = 0.6. Such lines delimit a line of continuous 
transitions that does not depend on a. The discontinuous lines 
end into a critical point. The plot shows that the stronger are 
the correlations between the activity of the nodes and their 
degree within each layer, the larger is the percolating phase, 
but at the same time the transition becomes discontinuous. 


In particular, we assume for simplicity that 

PB.(fc) = (33) 


and we take 


c{B) = coB^ 


(34) 


where co,a are two parameters determining the correla¬ 
tions between the degrees of the node and its activity. 
The particular choice of the functional form of Pb^ (k) in 
Eqs. (33)-(34) is dictated by the intention to model posi¬ 


tive correlations between the activity of the degree of the 
layers that have been observed in real dataset [7]. Real 
multiplex network analysis are nevertheless not sufficient 
to suggest the exact form of the correlations observed. 
Therefore this ensemble of correlated multiplex networks 
with heterogeneous activity of the nodes has been cho¬ 
sen in such way describe positive correlations between 
activities and degree in single layers, while keeping the 
model sufficiently simple to allow a number of analytical 
calculations. 

Averaging the message passing equation over this en¬ 
semble and assuming Na = N, Va, and B M, we have 
only one average message determined by the equation 


S = ^=pY. ' (l - e-=(^)'^)(35) 


B-1 


(B) 


Setting p = cop, the equation for x = cqct reads 


h(x) = 0 


(36) 


FIG. 7: Phase diagram of a multiplex network in the case 
of an activity distribution P{B) oc B~^ for B £ [2,100] and 
degree within each layer correlated with the activity of the 
node according to Eqs. (|33[) — (34l. From top to bottom. 


the dashed lines indicate discontinuous phase transitions for 
a = 0, a — 0.2, a = 1 and a = 1.5, respectively. Differently 
from the previous case, we do not have continuous transitions. 



X 

FIG. 8: Solutions of Eq. ( |37[ ) at the critical point for Bmin = 1 
and a — 0.2. 



FIG. 9: Position of critical points for a multiplex network 
in the case of an activity distribution P{B) oc B~^ for B £ 
]1,100] and degree on each la yer cor rela ted with the activity 
of the node according to Eqs.(331 — ([34|). 




















with 


h{x) = x-pJ2 ' (l - ■ 

This equation can be studied as a function of a, p and 
the parameters determining the P{B) distribution. In 
order to find the discontinuous phase transition, we set 
h{x*) = 0, h'{x*) = 0. The line of these discontinuous 
phase transition eventually stops at a critical point Xc, 
that can be calculated by setting 

h(xc) = h'{xc) = h"{xc) = 0. (37) 


The continuous phase transition can be found, instead, 
by imposing /i(0) = h'{Q) = 0. 

Let us now consider a power-law distribution of activ¬ 
ities P{B) oc B~^ for B G 100] and correlated 

degrees in every layer according to the degree distribu¬ 
tion given by Eq. (331. In Figure]^ and Figure]^ we show 
the percolation transitions as a function of the param¬ 
eter a for Bmin = 1 and Bmin = 2, respectively. As 
a general remark, for Bmin = 2, the transition appears 
to be always discontinuous, while for Bmin = 1 we have 
both continuous and discontinuous transitions, with a re¬ 
gion of coexistence between two percolating phases, and 
a critical point at the end of the line of discontinuous 
transitions. 

First, let us observe the case of Bmin = 1 (Fig. |^. 
We do not plot the line o = 0 as it reduces to the non- 
correlated case shown in Fig. where the phase transi¬ 
tion is always continuous. For a > 0, instead, it emerges 
a line of discontinuous phase transitions ending in a crit¬ 
ical point, which is determined by the Eqs. (37). The 
line of continuous phase transitions encounters the line 
of discontinuous transitions at a critical end point de¬ 
termined by the simultaneous occurrence of two minima 
of function h{x) (Eq. 37). Fig. shows the solution of 
equation (|3^ at the critical point. This phase diagram, 
therefore, is characterized by the presence of two perco¬ 
lating phases, separated by a short line of discontinuous 
transitions (Fig. [^. This coexistence region shifts to¬ 
wards larger values of (5 as a increases, but the critical 
point never joins the continuous line at a tricritical point 
at finite a. Figure]^ shows that Xc —t 0, without joining 
the continuous solution at a; = 0, for a —>■ oo. 

In the case Bmin = 2 (Figure [^, we observe a discon¬ 
tinuous phase transition for all values of a. This quali¬ 
tative difference is due to the peculiar role of the nodes 
which are active on a single layer [B = 1). In our model, 
this type of nodes do not need support from the other lay¬ 
ers (as they appear only in one of them), and therefore 


they drive a classical continuous percolation transition. 
In both cases, as a increases, the percolating threshold 
becomes smaller, as the intensity of the correlations be¬ 
tween the activities of the nodes and the degree of the 
nodes within each layer increases. Therefore the multi¬ 
plex networks is more robust because it can still have 
a mutually connected component even if the damage is 
significant. On the other hand, though, the phase tran¬ 
sition becomes discontinuous, and therefore the collapse 
becomes more unpredictable. 


V. CONCLUSIONS 


In this paper, we have characterized the robustness 
properties of multiplex networks in a new model that en¬ 
capsulates heterogeneous activities of the nodes, i.e. the 
possibility that each node is present only on a small frac¬ 
tion of layers in a multiplex, as seen in real world cases 
[7j. In this model, we employ a notion of mutual perco¬ 
lation where nodes must belong to a mutually connected 
component only on the layers where they are active and 
develop an analytical approach to calculate the size of 
the mutually connected giant component as a function of 
node damage and other parameters. We show that mul¬ 
tiplex networks with very broad activity distributions are 
more fragile than networks with more homogeneous dis¬ 
tribution of activities. 

We also investigate the role of correlations between 
the activities of the nodes and their degrees. We show 
that these correlations generally improve the stability of 
the percolating phase,and the multiplex network has a 
smaller percolation threshold, so the multiplex network 
becomes more robust. However, correlations also change 
the order of the phase transition, that becomes discon¬ 
tinuous. This provides an example on how correlations 
can typically reduce the fragility of multiplex networks, 
but at the same time they can make the system more 
unpredictable, as the transition becomes discontinuous. 
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